Bifurcations and Competing Coherent Structures in the 
Cubic-Quintic Ginzburg-Landau Equation I: Plane 

Wave (CW) Solutions 

Stefan C. Mancas 
and 

S. Roy Choudhury 
Department of Mathematics 
Embry Riddle Aeronautical University 
Daytona Beach, FL 32114 
mancass@erau.edu 
choudhur@longwood.cs.ucf.edu 

Abstract 

Singularity Theory is used to comprehensively investigate the bifurcations of the 
steady-states of the traveling wave ODEs of the cubic-quintic Ginzburg-Landau equa- 
tion (CGLE). These correspond to plane waves of the PDE. In addition to the most 
general situation, we also derive the degeneracy conditions on the eight coefficients of 
the CGLE under which the equation for the steady states assumes each of the possible 
quartic (the quartic fold and an unnamed form), cubic (the pitchfork and the winged 
cusp), and quadratic (four possible cases) normal forms for singularities of codimension 
up to three. Since the actual governing equations are employed, all results are globally 
valid, and not just of local applicability. In each case, the recognition problem for the 
unfolded singularity is treated. The transition varieties, i.e. the hysteresis, isola, and 
double limit curves are presented for each normal form. For both the most general case, 
as well as for various combinations of coefficients relevant to the particular cases, the 
bifurcations curves are mapped out in the various regions of parameter space delimited 
by these varieties. The multiplicities and interactions of the plane wave solutions are 
then comprehensively deduced from the bifurcation plots in each regime, and include 
features such as regimes of hysteresis among co-existing states, domains featuring more 
than one interval of hysteresis, and isola behavior featuring dynamics unrelated to the 
primary solution branch in limited ranges of parameter space. 
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1 Introduction 



The cubic complex Ginzburg-Landau equation (CGLE) is the canonical equation governing 
the weakly nonlinear behavior of dissipative systems in a wide variety of disciplines [1]. 
In fluid mechanics, it is also often referred to as the Newell- Whitehead equation after the 
authors who derived it in the context of Benard convection [H [2]. 

As such, it is also one of the most widely studied nonlinear equations. Many basic prop- 
erties of the equation and its solutions are reviewed in [21 H], together with applications 
to a vast variety of phenomena including nonlinear waves, second-order phase transitions, 
superconductivity, superfluidity, Bose-Einstein condensation, liquid crystals and string the- 
ory. The numerical studies by Brusch et al. O |6] which primarily consider periodic traveling 
wave solutions of the cubic CGLE, together with secondary pitchfork bifurcations and period 
doubling cascades into disordered turbulent regimes, also give comprehensive summaries of 
other work on this system. Early numerical studies |7] and theoretical investigations |8] of 
periodic solutions and secondary bifurcations are also of general interest for our work here. 

Certain situations or phenomena, such as where the cubic nonlinear term is close to zero, 
may require the inclusion of higher-order nonlinearities leading to the so-called cubic-quintic 
CGLE. This has proved to be a rich system with very diverse solution behaviors. In particu- 
lar, a relatively early and influential review by van Saarloos and Hohenberg [9], also recently 
extended to two coupled cubic CGL equations [ini E] , considered phase-plane counting ar- 
guments for traveling wave coherent structures, some analytic and perturbative solutions, 
limited comparisons to numerics, and so-called 'linear marginal stability analysis' to select 
the phase speed of the traveling waves. Among the multitude of other papers, we shall only 
refer to two sets of studies which will directly pertain to the work in this paper. The first 
of this series of papers [12]- [IS] using dynamical systems techniques to prove that the cubic- 
quintic CGLE admits periodic and quasi-periodic traveling wave solutions. The second class 
of papers [T7t [T8] , primarily involving numerical simulations of the full cubic-quintic CGL 
PDE in the context of Nonlinear Optics, revealed various branches of plane wave solutions 
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which are referred to as continuous wave (CW) solutions in the Optics hterature. More 
importantly, these latter studies also found various spatially confined coherent structures 
of the PDE, with envelopes which exhibit complicated temporal dynamics. In [17], these 
various structures are categorized as plain pulses (periodic), pulsating solitary waves and 
so on depending on the temporal behavior of the envelopes. In addition, the phase speed 
of the coherent structures may be zero, constant, or periodic (since it is an eigenvalue, the 
phase speed may be in principle also be quasiperiodic or chaotic, although no such cases 
appear to have been reported). Secondary complete period doubling cascades leading as 
usual to regimes of chaos are also found. This last feature for numerical solutions of the full 
cubic-quintic PDE is strongly reminiscent of the period doubling cascades found in |5l |6] for 
period solutions of the traveling wave reduced ODEs for the cubic CGLE. 

Motivated by the above, we begin a fresh look at the traveling wave solutions of the 
cubic-quintic CGLE in this paper. Besides attempting to understand the complex numerical 
coherent structures in [18], one other goal is to build a bridge between the dynamical systems 
approach in [12]- [IS] and the numerical one in [TTJ [18]. Given the importance of the cubic- 
quintic CGLE as a canonical pattern-forming system, this is clearly important in and of 
itself. However, a word of warning is in order here. Some of the features in [18] may well be 
inherently spatio-temporal in nature, so that a spatial traveling-wave reduction may not be 
sufficient to completely capture all aspects. Indeed, there is some evidence along these lines 

In this paper, we begin by using Singularity Theory |20j to comprehensively categorize 
the plane wave (CW) solutions which were partially considered numerically in [17]. In 
addition, we shall be able to identify co-existing CW solutions in all parameter regimes 
together with their stability. The resulting dynamic behaviors will include hysteresis among 
co-existing branches, as well as the existence of isolated solution branches (isolas) separated 
from the main solution branch. Subsequent papers will consider periodic traveling waves 
(traveling periodic wavetrains of the PDE), quasiperiodic solutions and homoclinic solutions 
(corresponding to pulse solutions of the PDE) and their bifurcations. 
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The remainder of this paper is organized as follows. Section 2 considers two formulations 
for the traveling-wave reduced ODEs for the cubic-quintic CGLE, as well as CW solutions. 
Section 3 quickly recapitulates the standard stability analysis for individual CW solutions. In 
Section 4, which is the heart of the paper, Singularity Theory is employed to comprehensively 
categorize all possible co-existing and competing plane wave solutions in general parameter 
regimes, as well as special cases corresponding to all possible quartic and cubic normal forms 
for singularities of codimension up to three. Section 5 considers the corresponding bifurcation 
diagrams as well as the resulting dynamical behaviors. 

2 Traveling Wave Reduced ODEs 
2.1 Reductions 

We shall consider the cubic-quintic CGLE in the form [9] 

dtA = eA + (6i + ic,)dlA - (63 - ic^)\A\^A - {b, - ic,)\A\^A (2.1) 

noting that any three of the coefficients (no two of which are in the same term) may be set 
to unity by appropriate scalings of time, space and A. 

For the most part, we shall employ the polar form used in earlier treatments O |9] of the 
traveling wave solutions of (12.11) . This takes the form of the ansatz 

A{x,t) =e-''^'A{x-vt) 

= e-''^'a{z)e''^^'^ (2.2) 

where 

z = X — vt (2.3) 

is the traveling wave variable and u and v are the frequency and translation speed (and are 
eigenvalues). Substitution of (12.21) / (12. 3p in (12.11) leads, after some simplification, to the three 
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mode dynamical system 



az = b (2.4a) 

2 (bie + Ciw)a + v{bib + Ciipa) - (6163 - 0103)0^ - (6165 - CiC5)a^ 
bz = aip -2- — 2 (2-4b) 



6? + c? 



_ _2# + cie + - - (&1C3 + &3Ci)q^ - (61C5 + b^ci)^ 

a ^ bl + c\ ^ > 

where ip = (pz- Note that we have put the equations into a form closer to that in [5J, rather 
than that in [9], so that (12.41) is a generalization of the traveling wave ODEs in [5] to include 
the quintic terms. 

For future reference, we also include the fourth-order ODE system one would obtain from 
(12. ip using the rectangular representation 

A{x,t) = e-'^^'Aix-vt) 

= e-''''[a{z)+tf3{z)] (2.5) 

with z given by (12. 3p . Using (12. 5p in (12. ip yields the system: 

-Ci5z + bi'jz = Ti (2.6a) 
biSz + ci7, = T2 (2.6b) 

where 7 = a', 5 = ' = d/dz, and ri/r2 are given below. This may be written as a first 
order system 

a' = 7 
P' = 6 

{bl + clh' = biTi + C1T2 

{bl + cl)6' = b,r2-c,r, (2.7) 
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with 



Ti = uP-v-f-ea + {ha + cg/?) (a' + /?') + (65a + cs/?) (a' + (2.8a; 



and 



= -ua ~v6-el3 + (63/? - Csa){a' + l3') + {^(3 - C5a)(a' + (5y. (2.8b) 

2.2 Fixed Points and Plane (Continuous) Wave Solutions 

From (12. 2p . a fixed point (ao,0,'?/'o) of (12. 4p corresponds to a plane wave solution 

A(a;,t) = aoe^('^o"-''*)+^^ (2.9) 

with 9 an arbitrary constant. 

The fixed points of (I2.4al) . fl2.4bp and (I2.4cp may be obtained by setting 6 = (from 



(I2.4ap ) in the right hand sides of the last two equations, solving the last one for if), and 
substituting this in the second yielding the quartic equation 

a^^x'^ + a^x^ + a2X^ + aix + = (2-10) 
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with 



X = a\ (2.11a) 
(&1C5 + h^c^f 

«4 = ^2 (2-llb) 

2(6iC3 + 63Ci)(6iC5 + 65C1) /o 11 ^ 

«3 = (2.11c) 

6^c? + 26163C1C3 - 265C?e 65t;2 + 6i(c| + 2c5w) + 201(650; -cge) 

63 2(610; - cie)(6iC3 + 63C1) /oil N 

tti = — + -^^ 2.11e 

61 O^f^ 

(cie-610;)^ e 
off"' 61 

Thus, with oq = a/x for x any of the four roots of (12.101) . we have a plane wave solution of 
the form ([2l]). 

The fixed points of the system (12.71) are given by 7 = 5 = and Fi = r2 = 0. They 
may be obtained by eliminating the a and /3 terms by solving Fi = r2 = simultaneously 
yielding: 

+ /^^ = 



or 



2 , n2 hu + c^e 



hc^ - 65 C3 

Resubstituting these into the Fi = r2 = yields only the trivial fixed point 

a = (3 = 0. (2.12) 

Thus, the system (12.71) has no non-trivial plane wave solutions. 

In the next section, we begin the consideration of the stability, co-existence and bifurca- 
tions of the plane wave states of (12. ip (the fixed points of (I2.4al) . (]2.4bp and (I2.4cl) ). 
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3 Stability Analysis for Individual Plane Wave 
Solutions 



In this section, we conduct a stability analysis of individual plane wave solutions using regular 
phase plane techniques. This was already done for the alternative formulation of the traveling 
wave ODEs given in [9]. We provide a brief derivation for our system fl2.4al) . (I2.4b|) and f l2.4cp 



for completeness and future use. However, a much more complex question is the issue of 
categorizing and elucidating the possible existence of, and transitions among, multiple plane 
wave states which may co-exist for the same parameter values in (12. ip (corresponding to the 
same operating conditions of the underlying system). Such behavior is well-documented in 
systems such as the Continuous Stirred Tank Reactor System [201 121] • For a system such 
as (12. ip and the associated ODEs (I2.4ap . ( ]2.4bp and (I2.4cp . the large number of parameters 



makes a comprehensive parametric study of co-existing states bewilderingly complex, if not 
actually impracticable. This more complex issue is addressed in the next section. 

For each of the four roots Xi,i = 1, ... ,4 of (I2.10p corresponding to a fixed point of 
(I2.4ap . ( I2.4bp and (I2.4cp or a plane wave ^/xi e*('^»^~'^*)+*^> , the stability may be determined 



using regular phase-plane analysis. The characteristic polynomial of the Jacobian matrix of 
a fixed point Xi = of of (I2.4ap . fl2.4bp and (I2.4cp may be expressed as 



A' + 5i\^ + + 53 = (3.1) 
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where 



5, = ^ (3.2a) 



^2 = ,2 ! 2 [3a^(ciC3 - 6163) + 5a'^(ciC5 - 6165) 
"1 + Ci 

+ (6ie + ciw) + 3{bl + cDij^ + v{v - 3ci^)] (3.2b) 
_ 4a2V^[(6iC3 + 63C1) + 2a2(6iC5 + ha)] 

+ (^2 ^^2)2 {^1^1^^^ - ^ a^(263(&? + C?) + 61(6163 - C1C3)) 

+ 0^(465(6? + cj) + h{b,b, - C1C5)) - 61 ((6ie + cico) - ij\hl + c?))] }, (3.2c) 

where the fixed point values (a^, ipi) = {^/Xi, ipi) are to be substituted in terms of the system 
parameters (§2). Note that ipi is obtained by setting and 6 = in the right 

side of (^m . 

Using the Routh-Hurwitz conditions, the corresponding fixed point is stable for 



5i > 0, 53 > 0, 5i52 - 53 > 0. (3.3) 

Equation (13.31) is thus the condition for stability of the plane wave corresponding to Xj. 

On the contrary, one may have the onset of instability of the plane wave solution occurring 
in one of two ways. In the first, one root of (13. ip (or one eigenvalue of the Jacobian) becomes 
non-hyperbolic by going through zero for 

53 = 0. (3.4) 

Equation (13. 4p is thus the condition for the onset of 'static' instability of the plane wave. 
Whether this bifurcation is a pitchfork or transcritical one, and its subcritical or supercritical 
nature, may be readily determined by deriving an appropriate canonical system in the vicinity 
of (13. 4p using any of a variety of normal form or perturbation methods [21]-j23]. 
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One may also have the onset of dynamic instabihty ('flutter' in the language of Ap- 
plied Mechanics) when a pair of eigenvalues of the Jacobian become purely imaginary. The 
consequent Hopf bifurcation at 

6^62 -63 = (3.5) 



leads to the onset of periodic solutions of fl2.4ap . fl2.4bp and fl2.4cp (dynamic instability 



or 'flutter'). These periodic solutions for a{z) and iplz), which may be stable or unstable 
depending on the super- or subcritical nature of the bifurcation, correspond via (12. 2p to 
solutions 

A{x,t) =a{z)e<^'^'''-''') (3.6) 

of the CGLE (12. ip which are, in general, quasiperiodic wavetrain solutions. This is because 
the period of ip and u are typically incommensurate. Eq. (13. 6p is periodic ii u = 0. We shall 
consider these wavetrains, including the derivation of normal forms, more general versions 
of the Hopf bifurcation, and stability, in a companion paper. 

Here, we change gears to address the more difficult question of the possible coexistence 
of, and transitions among, multiple plane wave states for the same parameter sets. 

4 Co-existing and Competing Plane Waves 

As mentioned earlier, for a multiparameter system like (12. ip . and the associated ODEs (I2.4ap . 
fl2.4bp and (I2.4cp . a comprehensive parametric study of co-existing states is forbiddingly 



complex, if not actually impracticable. Theoretical guidance is needed to determine all the 
multiplicity features in various parameter domains, as well as the stability of, and mutual 
transitions among, coexisting plane waves in each domain. 

In this section, we use Singularity Theory [20j to comprehensively analyze such multi- 
plicity features for (12. ip / (I2.4ap . fl2.4bp . (I2.4cp . In particular, we shall derive the existence 
conditions on the eight coefficients of the CGLE under which the steady state equation (I2.10p 
assumes either a. all possible quartic normal forms (the quartic fold, and an unnamed form). 
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or b. all distinct cubic normal forms (the pitchfork or the winged cusp) for singularities of 
codimension up to three. In addition, given that the most degenerate singularities or bi- 
furcations tend to be the primary organizing centers for the dynamics, we also consider the 
even higher codimension singularities leading to various quadratic normal forms. Clearly, the 
most degenerate singularities for a particular parameter set would 'organize' the dynamics 
in the sense that local behavior in its vicinity predicts actual quasi-global results. In fact, 
since we employ the actual governing equations, the ensuing results are not just locally valid, 
as is often the case, but they have global applicability. 
First, denoting fl2.10p as 



where g denotes the 'germ' and the are given in terms of system parameters by (12.111) . all 
points of bifurcation (where the Implicit Function Theorem fails) satisfy 



Given a germ satisfying (14. ip / (14. 2p . the general Classification Theorem in [20] provides a 
comprehensive list of all possible distinct normal forms to which it may be reduced for 
bifurcations of codimension less than or equal to three. 

For our g, which is already in polynomial form, it is particularly straightforward to reduce 
it to each of these normal forms in turn and this is what we shall do next. Following this, we 
shall consider the general form (14. ip itself. We start first with the possible distinct quartic 
normal forms viz, the 'Quartic Fold' and an unnamed form, and then proceed systematically 
to lower order normal forms. In the standard manner, the so-called 'Recognition Problem' 
or identification of each normal form yields certain defining conditions and non-degeneracy 
conditions and we check these first for each form. Each normal form has a well-known 
'universal unfolding' or canonical form under any possible perturbation [20]. This is so 




(4.1) 



0. 



(4.2) 
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under certain other non- degeneracy conditions (the conditions for the solution of the so-called 
recognition problem) which we next satisfy. Once the universal unfolding is established, we 
next need to consider the various parameter regions (for the parameters in the unfolding) 
where distinct behaviors for the solutions x occur. The boundaries of these regions are the 
so-called 'transition varieties' across which these behaviors change or are non-persistent. We 
consider these next. The final step involves detailing in each region delimited by two adjacent 
'transition variety' curves the bifurcation diagram for x, i.e., the possible co-existing steady 
states of f l2.4ap . ( ]2.4bl) and f l2.4cp (or plane waves of (12. ip ) and their stability. 



4.1 The Quartic Fold 

We perform the steps mentioned above for the first quartic normal form, viz. the Quartic 
Fold 

hi{x,X) = ex^ + 6X. (4.3) 
Clearly, our germ (4.1) has this form for 

a4 = e, 

as = 0(2 = C(i = 

ao = S\. (4.4) 
For the normal form (14.31) . the universal unfolding is 

Gi (x. A) = ex^ + 5X + ax + (3x^ (4.5) 

with defining conditions 

Qxx = Qxxx = 0, (4.6) 
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non- degeneracy conditions 

and provided the condition for the solution of the recognition problem 



9x 


9\x 


9\xx 






Gia 


G\ax 


Glctxx 




(4.7) 


Gli3 


Gipx 


Gipxx 







is satisfied. Given (14.1 p and (14. 4p . the conditions (14. 5 p are automatically satisfied, while (14. 7p 
yields the condition 

5^0. (4.8) 
The transition varieties across which the (A,a;) bifurcation diagrams change character 

are: 

i. The Bifurcation Variety 

= {aeR^ : (x, A) such that G = Gx = Gx = Q at (x, A, a)}. (4.9) 

ii. The Hysteresis Variety 

'K= {deR^ : {x, A) such that G = = = at (x, A, a)}. (4.10) 

and, 

iii. The Double Limit Variety 

V = {a e R^ : {xi,X2,\),xi X2 such that G = Gx = Q at (xj. A, a), z = 1, 2}. (4.11) 
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We compute these here since the derivations are not given in |20|. For S, we need 

Gi^ = 4ex^ + a + 2/3x = 

and 

Gu = 5 = 0. 

However, (5 7^ by ( I4.8p . and hence the bifurcation set is just the null set 

'3 = 0. (4.12) 

For "K, we need 

Gi^ = 4ex^ + a + 2/3x = 
Gi^^ = VlZx^ + 2/3 = 



which yield 

Similarly, using f l4.1ip . it is straightforward to derive the double limit set 



^=<{^y = -(iX, /3<0^ (4.13) 



2) = {a = 0,/3 < 0}. (4.14) 

In the (a, /3) plane, the (A, x) bifurcation diagrams change character across the curves f l4.12p . 
f l4.13p and fl4.14p . so that there are different multiplicities of steady-states in the regions they 
de-limit. We shall consider this in detail in the next section. 
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4.2 A Second Quartic Normal Form 

Repeating the above steps for the other possible distinct normal form 



h2{x, X) = ex + 6Xx 



(4.15) 



our germ (14.11) takes this form for 



0^4 = e 

ttS = tt2 = tto = 
tti = 6X. 

For the normal form (I4.15p . the universal unfolding is 

G2{x, X) = ex'^ + SXx + a + (3X + 

with defining conditions 

9xx 9xxx 9\ 0, 



(4.16) 



(4.17) 



(4.18) 



non-degeneracy conditions which are automatically satisfied, and the solution of the recog 
nition problem yielding the condition 









9x\ 





9xxxx 





9xx 


9\x 


9Xxx 


9Xxxx 


G2a 


G2ax 


G2aX 


G2axx 


G2axxx 


G213 


G2j3x 


G2(5X 


G213XX 


G213XXX 


G2-1 


G2-yx 


G2-yX 


G2'yxx 


G2'yxxx 



7^0. 



(4.19) 



For CT /f lCTD . Km is satisfied, while ( ICTll yields 

eS 7^ 0, or aia^ 7^ 0. 
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(4.20) 



We derive the transition varieties for this case since derivations are not provided in [20^ 

For 



G2x = 4ex^ + 6X + 27X = 
G2x = Sx + /3 = 



which, together with (I4.17p . yield 



a3:^ + ^ + a = 0. (4.21) 



For JC: 



and 



G2XX = 0^7 = —6ex 



Together, these yield 

= -87V275^e. 
Using these in fl4.17p yields the hysteresis curve: 



, 72 \^ 87^/32 



Similarly, the double limit curve 2) is: 

2? : 4a = 7^, 7 < 0. (4.23) 

In the next two subsections, we summarize similar results for the two distinct cubic 
normal forms, but omit the details. Then we briefly mention the four possible quadratic 
normal forms for even more degenerate cases, before concluding the section with the general, 
least degenerate case. 
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4.3 The Pitchfork 

For our germ (14. ip to have the cubic normal form for the well-known pitchfork bifurcation 



/i3(x,A) = ex^ + (5Ax (4.24) 



we require 



0^4 = ^2 = ao = 

as = e, ai = 5\. (4.25) 



This will have a universal unfolding [20] 



Gg = tx^ + 5Xx + a + I3x^ (4.26) 



provided 











h^xxx 









h^Xxx 




G^ctx 


GsaX 


G^axx 


G3I3 


GsiBx 


Gspx 


GsfSxx 



= Sy^O. 



The well-known transition varieties, generalized to our notation, are: 

!B : a = (4.27) 

:K : a = ^^jnt (4.28) 

2) : 0. (4.29) 
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4.4 The Winged Cusp 

The other distinct cubic normal form 

h^ix, A) = ex^ + (5A^ 

requires 

0^4 = ^2 = = 

as = e, ao = 5A^. 

This has a universal unfolding [20] 

G4(x, A) = ex^ + dX^ + a + l3x + 'yXx 

provided 









hix\ 


^Axxx 





hiXx 




h^Xxx 




G^ax 


GiaX 


G^axx 


G413 


Gai3x 


Gai3X 


Ga/^xx 



= -12Se ^ 0. 



The transition varieties, for our unfolding G4, are: 

S : a = 2x^ - /3 = -3x^ + -f^x/2 

: aj^ + = 0, a<0 
V : 0. 
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4.5 Quadratic Normal Forms 

Since our system of ODEs has many parameters, we may clearly have more degenerate 
(higher codimension) cases corresponding to any of the distinct quadratic normal forms 

h5{x,X) = ex^ + 5\ (4.36) 

/ig(a;,A) =e(x2 - A^) (4.37) 

/i7(x, A) = e(x2 + A^) (4.38) 

h8ix,\) =ex^ + (4.39) 

or 

hg{x,X) = ex'^ + 6X^ (4.40) 

Each of these is obtained by matching our germ (14. ip to the appropriate form, with the 
defining and non-degeneracy conditions automatically being satisfied (because (14.11) is poly- 
nomial). Solving the recognition problem [20], the corresponding unfoldings are respectively 

G5{x,X) = ex'^ + 6X (4.41) 
Gej{x,X) = e{x'^ + 6X'^ + a) (4.42) 

(with 5 < for (lOTll and 5 > for (ICTl) ) 

G8ix,X) = ex'^ + SX^ + a + /3X (4.43) 
Gg{x,X) = ex'^ + 6X^ + a + (3X + -fX'^ (4.44) 

with determinant conditions [20] for the cases f l4.43p and (14.440 which may be straightfor- 
wardly enforced as in previous cases. The S, "K, and D curves for these cases are straight- 
forward generalizations of those given in [20], and they may be derived as for the quartic 
and cubic cases. 
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4.6 General Case 

Finally, we include the most general possibility where, for arbitrary parameters in the CGLE 
fl2.ip . we have the germ f l4.ip with all non-zero. Treating f l4.ip itself as the unfolding, 
with aQ the bifurcation parameter A, the transition varieties in the (01,02) plane are: 



S : (4.45) 

ai = 8a4X^ + Saax^ (4.46) 
D : identical to "K (see Theorem 1 below) (4.47) 



Theorem 1. The Double Limit Variety for ( 14. ip is identical to the Hysteresis Variety of 
Proof. Using (14. ip and (14. lip . D is defined by the equations 



G'(xi,A) = (4.48a) 

G'(x2,A) = (4.48b) 

G,(xi,A) = (4.48c) 

G^(x2,A) = 0. (4.48d) 



Canceling the trivial solution xi = X2, the equations obtained from the difference of (I4.48ap 
and f l4.48bp . and of (I4.48cp and f!4.48dp . yield respectively 

ai = —a{2b — a^)a4 — has — aa2 (4.49a) 
^2 = -^(46a4 + Saas) (4.49b) 
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where a = Xi + x^, b = xl + X1X2 + x\. Using (14.49bp in fl4.49ap . these yield 



ai = ^-^ hja^ + aai (4.50a) 

a2 = —-^Cis — 2ba4. (4.50b) 

The equations (ICTal) and f l4.48bD may be considered to define the bifurcation parameter 
ao which we do not require here. However, f l4.48cl) and ( ]4.48d[) independently define a\ and 
thus far only their difference has been used. In order to incorporate ai, we consider the sum 
of fl4.48cp and fl4.48dp written in terms of a and b as: 

4a4a(36 - 2a^) + 'ia^i^h - c?) + 2a2a + 2ai = 0. 

Using fl4.50p in this equation and simplifying yields 

h = — . (4.51) 
4 

Using this in f l4.50p yields the parametric equations for D: 

ai = a [— — h 0:40] (4.52a) 
a2 = — — {0:3 + 0:40) . (4.52b) 

The re-parametrization a = 2x puts this into exactly the form fl4.46p of the hysteresis 
variety, thus proving the claim. □ 

Note that the "K curve is parametrized in terms of x (with 03,04 being chosen values). 
Also, given the non-degenerate nature of this general case, it is not surprising that there is 
only one distinct transition variety. 
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5 Bifurcation Diagrams and Effects on the Dynamics 

Having mapped out the S, !K, and D curves for the various possible distinct quartic and cubic 
normal forms, we now proceed in this section to consider the various bifurcation diagrams 
in the regions which they define in {a, /3) space. These will then give us the multiplicities 
and stabilities of the various co-exisitng steady states of fl2.4ap . f l2.4b|) and f l2.4cp (or plane 



wave solutions of f l2.ip ) in each region. In turn, these also enable us to consider dynamic 
features of the plane wave solutions. The dynamics will include hysteretic behaviors among 
co-existing plane waves. We will also find regimes of isolated plane wave behavior, both for 
a plane wave branch which co-exists with other branches but cannot interact with them, as 
well as those which actually occur only in isolation. 

We first list examples of representative sets of parameters for which we may have the 
various degenerate cases considered in Section 4. 

a. For the Quartic Fold of Section Section 4.1, typical parameters are: 

i. bi = 0.0845, 63 = -0.0846, 65 = 0.0846, ci = C3 = -C5 = 1, e = 0.5, v = 0.1, 
cu = 0. 

ii. bi = h = 0.01696, 63 = -0.0206, Ci = 1, C3 = 1.25, C5 = -1, e = 0.5, v = 0.1, 
u = 0. 

b. For the quartic normal form of Section 4.2: 

i. bi = 2.035, 63 = 29.274, 65 = 9.8496, Ci = -0.1, C3 = -1, C5 = 0.08, e = 0.3, 

V = 0.3, u = 0.1. 

c For the Pitchfork case of Section 4.3: 

i. 61 = 0.0904, 63 = 0.0679, 65 = 0.1811, Ci = -0.4, C3 = 0.35, C5 = 0.8, e = 0.2, 

V = 0.01, u = -0.9. 

ii. bi = 0.0904, 63 = 0.0823, 65 = -0.1808, Ci = -0.4, C3 = 0.35, C5 = -0.8, e = 0.2, 

V = 0.01, u = -0.9. 
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Figure 1: Transition varieties for the winged cusp f l4.32p with e = 1 = 6 for the cases 7 < 0, 
7 = 0, and 7 > 0, respectively. IK is in sohd hues, and S is dashed. 

d. For the Winged Cusp of Section 4.4: 

i. 61 = 0.000923, 63 = +.00005548, 65 = 0.0013, Ci = 0.5, C3 = -0.03, C5 = -0.7, 
e = 0.01, V = 0.1, to = 0.15. 

For the winged cusp unfolding (14.321) in the particular form 

Gi{x, X) = + X'^ + a + /3x + -fXx = 0, 

the transition varieties fl4.33p and (14.341) are shown in the (a, (3) plane in Figure [T^,b,c 
for 7 < 0, 7 = 0, and 7 > 0, respectively. They divide the space into seven distinct 

regions. As mentioned earlier, the (A, x) bifurcation diagrams are isomorphous or 'persistent' 
or of similar form within each region, and they change form across the transition varieties (or 
'nonpersistence' curves) as one crosses into an adjacent region. The representative bifurcation 
diagrams in each of the seven regions are shown in Figure [2], and they give us a comprehensive 
picture of the co-existing plane wave solutions of (12. ip and their stability (given by the 
eigenvalues of the Jacobian, or here just the sign of Gx) in each region. Hence, as we shall 
consider next, one also has a clear picture of the ensuing dynamics from the plane wave 
interactions. 

First, note the mushroom shaped bifurcation diagram in Figure [2)d for region 2 of Figure 
[U Clearly, there are two distinct ranges of A (at the two ends of the mushroom) where three 
plane waves co-exist (with the central one being unstable). Thus the dynamics exhibits 
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Figure 2: The (A,x) bifurcation diagrams in the regions 1-7 of Figure [Tj respectively. 
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hysteresis. For instance, if A is decreased from large values, one stays on the lower branch 
until point A before jumping to the upper branch. If A is then increased, one stays on the 
upper branch until B and then jumps back down to the lower one. Similar hysteresis occurs 
in regions 4 through 7 of Figured] as seen in the corresponding bifurcation diagrams of Figure 
[2jl-g. In each case, hysteresis occurs between the upper and lower fixed points in the range 
of A with three co-existing solutions (the central one is always unstable). 

Another type of behavior is the isola, i.e., an isolated branch of solutions unconnected 
to the primary solution (the one at A — ±00). Such isola type behavior is seen in Figure 
[2l:;,e,g corresponding to regions 3, 5, and 7 of Figure [1] . In each case, the isola co-exists 
with the primary solution branch and is the chosen branch or not according to the initial 
conditions. However, once chosen, the dynamics is on the isola while A is in the domain of 
its existence once we leave this domain, the solution cannot jump to the primary branch and 
just disappears. 

Next, we consider the normal form (]4.15p in Section 4.2. Considering the unfolding (14.17^ 
in the particular form 

G2{x, A) = - Ax + a + /3A + = 0, 

the transition varieties fl4.2ip . fl4.22p . and fl4.23p are shown in Figure [3^,b,c for the cases 
7 > 0, 7 = 0, and 7 < respectively. Note in particular, a significant correction to [3] in 
the "K curve of Figure [3b. The "K curve (14.220 represents a pair of straight lines in the {a, (3) 
plane, rather than the incorrect form 

"+12i+275^ = ° 

in [20]. In Figure [3}:;, one consequence is two new regions or domains 13 and 14 of the {a, /3) 
space. Also, the bifurcation plots in the domains 3, 4, 5 and 8 are significantly modified 
from those given in [20] for the corresponding regions. 



25 



Figure 3: The transition varieties for the quartic normal form (14.171) with 5 = —1 for the 
cases 7 > 0, 7 = 0, and 7 < 0, respectively. "K is in solid lines, "B is dashed and the double 
limit curve 2) is in fine dashing. The regions 1-14 which they delimit are shown. 

The bifurcations plots in the fourteen regions in Figure [3t (Figure [3K,b feature only some 
of the regions) are shown in Figures H] and [51 Note that there are no regions of isola behavior. 
In regions 3, 4, 5, and 8, there is only one branch of solutions, rather than two as shown in 
Figure m (case 10) of [20]. Of these, the segments BC and DE are unstable in cases 3 and 
5, so that the hysteretic behavior of the solutions will consist of transitions from the stable 
plane waves on branch AB to those on branch CD as A is increased past point B, and a 
reverse transition when it is decreased through C. Similarly, in regions 4 and 8 where only 
segment BC is unstable, hysteresis occurs with a transition from the plane wave on branch 
DE to branch AB if A is decreased through D, a transition from branch CD to branch AB 
when A is decreased through D, and a transition from CD to either AB or DE (depending 
on system bias, noise et cetera) as A is increased through C. Analogous hysteresis behavior 
is clearly possible in regions 7 and 11, while regions 9, 10, and 12 feature hysteresis between 
co-existing stable plane wave solutions on distinct solution branches. In the two new regions 
13 and 14 of Figure [3b (which were missing in [20j), the bifurcation plots in Figures [Sjn 
and[5)i show only two co-existing plane wave solutions in each A range, unlike the adjacent 
regions 5 and 8 of Figure [3b where the bifurcation plots Figures [lb and [Di have A ranges 
with four coeval solutions. 

For the very degenerate cases discussed in Section 4.5 and corresponding to quadratic 
normal forms, the corresponding transition varieties as well as the bifurcation plots and 
resulting dynamics in the regions of {a, /3) which they delimit may be deduced from the 
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Figure 4: Bifurcation diagrams in the regions 1-8 of Figure Eb. 




Figure 5: Bifurcation diagrams in the regions 9-14 of Figure [31:;. 
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relevant cases in Figures 4.1-4.3 of [20]. In particular, cases f l4.42p . fl4.43p and fl4.44p show 
isola, hysteresis, and double isola behaviors respectively. 

Let us consider the general case in Section 4.6 next. Since the results are entirely new, 
we need to consider the issue of the stability of the {X,x) = {ao,x) bifurcation diagrams in 
various regions of the (ai, 0^2) plane. Using (14. ip and the Chain Rule, 

dG _ dG _ dG dx dG _ ^ 
dX dao dx dX dX 

so that 

^ = L_ (5 1) 

dx (dx/dX)' ^ ■ ^ 

Thus, the Jacobian and its eigenvalue Gx (these are identical for a one- dimensional system 
such as (14. ip ) are negative, and the corresponding fixed-point branch of the {X,x) plane is 
stable, for segments of the bifurcation plot where 

dx 

— > 0, stable. (5.2) 
dA 

Conversely, segments with ^ < are unstable. 

Finally, let us consider the dynamics and interactions of plane waves for the most gen- 
eral case of Section 4.6. The coincident transition varieties "K and D in (14.461) / (I4.47P are 
shown in Figures [6^-l6h for various combinations of (0:3,0:4) values. As is readily apparent, 
the configurations in Figures are the independent ones corresponding to centered and 
off-centered cusps and a parabolic variety curve respectively - the other cases are simple 
reflections of these. For Figure [6^ with (03,04) = (0, 1), the transition variety divides the 
(oi, 02) space into two distinct regions 1 and 2. The bifurcation plots in the two regions are 
shown in Figures [7^ and [7)d. As per (15. ip . the segment(s) with ^ > are stable, so that 
there is a unique stable plane wave for Figure [7^ in region 1 of Figure [6^. By contrast, there 
are co-existing stable plane wave states in regions BC and DE of Figure Uh (for region 2 of 
Figure [6^). Thus, hysteretic dynamics occurs with a transition from BC to DE as A = Oq is 
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Figure 6: Transition varieties for the general case (14. ip treated in Section 4.6. There is no 
S curve, and !K and 2) are coincident. The regions they dehmit are shown. The figures 
correspond respectively to (03,04) values (0,1), (1,1), (1,0), (—1,1), (—1,0), (—1,-1), 
(0,-1), and (1,-1). 

decreased through C, and a reverse transition as A is increased on DE through D. For Figure 
[6h with (03,04) = (1, 1), the bifurcation plots in regions 1 and 2 are shown in Figures |8^ 
and [8]d respectively. Once gain, per (15. ip . the segments of these plots with positive slope 
correspond to stable plane waves. Thus, only the segment corresponding to DE in Figure 
is a unique stable plane wave solution in region 1 of Figure (6)3. For region 2 of Figure |6]d. 
Figure [8]d shows hysteresis between the stable plane wave branches BC and DE. For regions 
1 and 2 of Figure Et corresponding to (03,04) = (1,0), the bifurcation plots in regions 1 
and 2 are shown in Figures |9^ and[9]D. For the former, as per (15. ip . no stable plane waves 
exist. For the latter, there is a unique stable plane wave solution in the range of A(oo) 
corresponding to segment BC. 

Finally, for the sake of completeness, we mention an alternative interpretation of the gen- 
eral case using Catastrophe Theory (see [2U] for a discussion of the connection between 
this and the Singularity Theory approach). Treating (14. ip in a manner analogous to the 
Cusp Catastrophe, 

Gx = 4o4X'^ -|- 3o3X^ -|- 202X -|- Oi 

= 4o4(x^ + Fax^ + Fix + Fq) = (5.3) 
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Figure 7: Bifurcation diagrams in regions 1 and 2 of Figure [6^, respectively. 
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the transition cusp curve between domains with one and three real solutions is given by 



+ = 0. 



(5.6) 



For (as, a^) = (1, 1) corresponding to Figures [6b and[8l the catastrophe surface (15. 3 p showing 
regions of one/three real solutions in the (ai, 0:2) plane shown in Figure [Hti'- Figure [Hb shows 
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A 

Figure 8: Bifurcation diagrams in regions 1 and 2 of Figure [6]d, respectively. 

the corresponding cusp surface (15. 6p in (ai,a2) space. As mentioned, [20] discusses the 
relationship between these plots and the Singularity Theory plots given in Figures and [S] 
for this case. 

In concluding, we have comprehensively analyzed the co-existing plane wave solutions 
in various parameter regimes for the CGLE (12.1 1) . This includes transitions among co- 
existing states involving up to two domains with hysteresis, isolated parameter regimes with 
isola behavior, and the resulting dynamics. We should also stress that, since our governing 
equation (14. ip is of polynomial form, all the results in Sections 4 and 5 are globally (and 
not just locally) valid in their respective regimes, as of course are the results for the general 
case. 
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